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Abstract 

We investigate the properties of the neutron star with relativistic mean field models. We incorpo¬ 
rate in the quantum hadrodynamics and in the quark-meson coupling models a possible reduction 
of meson masses in nuclear matter. The equation of state for neutron star matter is obtained and 
is employed in Oppenheimer-Volkov equation to extract the maximum mass of the stable neutron 
star. We find that the equation of state, the composition and the properties of the neutron stars 
are sensitive to the values of the meson masses in medium. 
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I. INTRODUCTION 


The state of matter at extremely hot or dense conditions is one of the most fundamental 
questions in physics. Terrestrial nuclei of heavy elements provide a glimpse at such extreme 
states, but there are far-more-extreme states of matter in our universe. A neutron star, which 
may be regarded as a huge nucleus, is one example of such an extreme state of matter. The 
state of matter can be characterized by the equation of state (EoS). From the microscopic 
point of view the EoS is determined by how the constituent particles of the matter interact 
with each other. Therefore the state of matter at extreme conditions can be understood 
when we have enough knowledge about the properties of the constituent particles and how 
they interact in such conditions. 

In the early 90’s, Brown and Rho proclaimed the scaling property of hadron masses 


Q- 


They 


in dense medium in terms of the scale invariance of the effective lagrangian 
showed that the ratios of the in-medium masses of the nucleon, a-, a;- and p-mesons to 
their masses in free space are approximately equal to each other at around the nuclear 
saturation density. Afterwards, changes in the hadron masses in medium were calculated 


in the framework of QCD-sum rule 


Refs 


. 0] and P|, 


□ 


In 


and quark-meson-coupling (QMC) models 
masses of p and u mesons were calculated at the nuclear saturation density, 
and similar amount of mass reduction from the free mass, 'm* pu] l m p,oj — 0.8 was predicted. 
Experimentally, microscopic properties of hadrons as well as the states of matter at hot or 
dense environment can be probed in relativistic heavy-ion collisions, and possible meson mass 
reductions in nuclear matter were reported through the dileptonic decay of p and u mesons 


in the CERES/NA45 and KEK-PS E325 |5[ experiments. Also, recent astronomical 
observations of binary systems yield the mass of compact objects as Mj 0 75 i+i 807 = (2.2 ± 


O.2)M 0 jg| and M 4U i 7 oo -37 = (2.44 ± 0.27)M o |7|, of which the former is believed to be a 
neutron star. Such a large mass of the neutron star requires a substantially stiff EoS. 

Motivated by these theoretical and experimental studies, in this work we take into ac¬ 
count the reduction of meson masses in the nuclear models such as quantum hadrodynamics 
(QHD) and QMC, and apply these models to the neutron star matter. The nuclear models 
are calibrated to the bulk properties of nuclear matter at the saturation density, such as 
saturation density p 0 (= 0.17 fm 3 ), binding energy E B (= 16.0 MeV), symmetry energy a sym 
(= 32.5 MeV) and compression modulus K (200 ~ 300 MeV). We take into account the 
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effect of meson mass changes in the following ways. First, we use the meson masses in free 


space for the nuclear models considered 
adopt the Brown-Rho (BR) scaling law 


lere as is usually done in most works. Secondly, we 


a, 


m 


N 


m N 


rri„ 


m n 


m„ 


m u 


m; 


nir 


( 1 ) 


and assume a density-dependent scaling function for the meson masses p]. The nuclear 
models that employ the scaling function are denoted by the abbreviation “S”. Thirdly, 
in using the QMC type models we may treat the heavy mesons (p and u) as meson bags 
(denoted by “MB”) composed of a quark and an anti-quark jfj]. In such cases, heavy mesons 
are treated in consistent with the description of the nucleon in the framework of the QMC. 
After calibrating different nuclear models for nuclear matter we apply them to the neutron 
star matter and see how the reduction of meson masses in matter may influence the EoS of 
the neutron star matter and the properties of the neutron star. 

In Sec. 2, we briefly describe the nuclear models considered in this work and fix the 
parameters of the models and show the resultant properties of the nuclear matter at the 
saturation density. In Sec. 3, we apply the various nuclear models fixed in Sec. 2 to the 
neutron star matter and discuss the results and implications for the properties of the neutron 
star. A summary is given in Sec. 4. 


II. MODELS FOR NUCLEAR MATTER 

In this work we will consider 5 different models. We refer to them as QHD, QHD-S; 
MQMC, MQMC-S and MQMC-MB. All the models have three coupling constants g a N , 
g U N, and g p ^ (or g q , g q , and g q p ) for the interaction of nucleons (or quarks) with cr-, to- and 
p-mesons, respectively, a- and uncoupling constants are fitted to reproduce the given values 
of p 0 and E B , and the coupling to p-meson is determined by a sym . Details of each model are 
explained in the following subsections. 


A. QHD 

The lagrangian for QHD in the mean held approximation reads 

l qhd = $N[i'fd ll - (m N - g aN a) - g^ N 7 0 u 0 ~ 7; 9 p n 7° ^03 t 3 ]^ n 
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~\ m l 0-2 _ \ m N b (dcrN O ') 3 - (fjtjN CT ) 4 

+\ m l^o+ \ m lbl 3 . (2) 

The masses uin, and m p refer to the values in free space, 939, 550, 783 and 770 

MeV, respectively. The effective mass of the nucleon is defined as 

m* N ( QHD) = m N - g aN a. (3) 


The meson helds are determined from the equations of motion in the mean field approxima¬ 
tion ; 
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E Imk% = (Pp-Pn), 


N=n,p 

where Rv 3 is 1/2 (-1/2) for the proton (neutron). 

1 * ^ pk N m 


2 m 2 


Ps = 
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k dk 


(4) 

(5) 

( 6 ) 

(7) 


71-2 N=n,p ® \Jk 2 + m*N 

is the scalar density and k n is the Fermi momentum of the nucleon at a given density. 
Cubic and quartic self-interaction terms of the cr-meson are included in Eq. © to produce 
reasonable K and m* N values at the saturation density. The coupling constants fixed to 
produce the saturation properties and the resulting m* N and K are listed in Tab. |H 


B. QHD-S 


We now consider the meson mass changes in matter and adopt the model proposed in 
Ref. B where the BR-scaling law is incorporated in the original QHD. The model lagrangian 
reads 


1 


Zqhd-s = - ( M n - 9oN cr) ~ g* N 7 - ^9 p n 7 b 03 r 3 ]^ N 

~\ m *a + \mf bis- 


The BR-scaling law is parametrized as 
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Model 

9aN 

9ujN 

dpN 

b x 100 

c x 100 

y 


m* N /m N 


K 

QHD 

8.11 

8.36 

7.85 

0.3478 

1.328 



0.773 

1.0 

310.8 

QHD-S 

5.30 

15.30 

7.52 



0.28 

0.31 

0.667 

0.781 

264.5 


TABLE I: Parameters and the saturation properties for the QHD and QHD-S models. The sub¬ 
script M for m* M /rriM refers to a, u and p mesons. 

The effective mass of the nucleon is defined as 


m w(QHD-S) = — g„ N a. 


( 10 ) 


(Note the use of M^ in place of of Eq. ©.) Coupling constants may change with nuclear 


matter density. Such density c 


considered in previous 


ependences of coupling constants as well as masses have been 


111 IT 121 and very recent work o : 

relativistic mean field theory. In Ref. jgj] , g^N is assumed to vary with density to satisfy 


9un/9uN — at around the saturation density. g^ N scaling is expressed as 

** = [i + Z D) 1 . (11) 

9ujN V Po / 

The parameters for the QHD-S model and the resulting properties at the saturation are 
summarized in Tab. |H One can see that both QHD and QHD-S models give K values in 
the acceptable range. rrjX^/mM turns out to be 0.781, which is close to the values obtained 
by previous studies |2j, 3j. rR N /rri N from QHD-S is smaller than that from QHD because of 
the scaling in Eq. ©• 


C. MQMC 

In the QMC model [15|, interactions between the nucleons (or bags) are mediated by 
the exchange of mesons which couple with the quarks in the nucleon bags. The mean held 
lagrangian for non-overlapping spherical bags in dense matter may be written as 


-MF 

7 QMC 


= ^ q [iYdn ~ {m° - g q a a)- gl 7 0 u 0 - -g q 7 0 b 0 3 r 3 - B] x 9 V (R - r)ijj q 


^ 2 - 2,1 2 - 2,1 2 - 2 
+ 2 mLuUJ ° + 2 mpb ° 3 ’ 


( 12 ) 


where g q M (q = u, d; M = a, ui, p) is the coupling constant for the quark-meson interaction. 
B is the bag constant and R is the bag radius within which quarks are conhned. The effective 
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mass of the nucleon in the QMC is given by 


m^(QMC) = 


rpN 
^bag 


Zj\r 4 


= y al - — + -nR 3 B, 

4" R R 3 


!J S = (m' q - m“ - j». 


(13) 

(14) 

(15) 


x q is the eigen energy of the quarks in the bag determined by the boundary conditions at 
r = R. Zn is a phenomenological constant that incorporates the effects not explicitly taken 
into account, including zero-point motion. 

The scalar and the vector potentials obtained from the QMC model were much smaller 
than those from the QHD, which would result in too weak a spin-orbit potential to explain 
the spin-orbit splittings in finite nuclei and spin observables in nucleon-nucleus scattering. To 
circumvent these shortcomings, density dependent bag constants were introduced [3]. The 
QMC model with a density dependent bag constant is called the modified QMC (MQMC) 
model. We employ the direct coupling model of Ref. [lil where the bag constant reads 


B = B 0 (l-gZ-= 


A <J V 


8 rriN 


> 


( 16 ) 


a value is determined from the self-consistency condition (SCC) 


a = 3 




■Ps 


CnA ) + 


B_ / _ igja 
g% m* N 9 iriN \ 8 tun _ 
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with 


t?N 

Cn{o) = 
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N 


n n 


E? ag R, 


S(a) + 


ml 


E 


N 


and 


S{&) = 


Qq/2 + Rm* (fi, - 1) 


Qq (Q q - 1) + Rm*/ 2' 

B 0 and Zjy are fitted to reproduce the free nucleon mass with the stability condition 


dm* N 


dR 


= 0 


(17) 


(18) 


(19) 


( 20 ) 


R=Rq 


where m* N is evaluated by Eq. (fTT)l . Choosing R 0 = 0.6 fin, we obtain B^ 1 = 188.1 MeV 
and Z N = 2.030. g® and 8, g q and g q are adjusted to produce saturation properties and 
reasonable values of K and m* N . g q p is fitted for the symmetry energy to be reproduced. 
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D. MQMC-S 


As a simple way to incorporate the reduction of meson mass in medium, we assume 



( 21 ) 


The scaling parameter y is determined to satisfy Eq. CD at around the saturation density 


with m* N given by Eq. (THfl) . The lagrangian for this model can be easily obtained by replacing 
the meson masses in Eq. ill 21 ) with the medium-modified values in Eq. ( 12111 . Since the SCC 
is obtained from the minimization of the energy density with respect to the variation of a, 
the SCC of the MQMC-S is of the same form as that of MQMC, but with m a replaced by 
m* given by Eq. m■ The model parameters fixed to produce saturation properties are 
listed in Tab. HU 

E. MQMC-MB 

In the framework of MQMC, we may treat vector mesons as MIT bags composed of a 
quark and an anti-quark. We refer to such a model as MQMC-MB. The parameters Bq, 
Ro and Zm (M = p, u) for the meson bags can differ depending on the mesons, but for 
simplicity we fix Bo and Ro as those values for the nucleon and treat only Zm as a parameter 
for each meson 191. The effective mass of each meson can be written as 



( 22 ) 


(23) 


Z M is fixed so that m* M defined as in Eq. (12211 recovers the mass of the u- and p -mesons in 
free space. Z u = 0.7904 and Z p = 0.8154 are obtained. The lagrangian for the model can be 
written by replacing the meson masses in Eq. (11211 with m* M in Eq. 1221 ) . In the MQMC-MB 
meson masses are functions of a, which modifies the SCC as 



7 




















Model 

9$ 

9l 

9a 

9 q P 

y 

m* N /m N 


K 

MQMC 

1.0 

2.71 

6.81 

7.89 


0.783 

1.0 

285.5 

MQMC-S 

1.0 

2.31 

5.51 

6.09 

0.28 

0.758 

0.781 

591.5 

MQMC-MB 

1.0 

1.77 

5.44 

8.15 


0.852 

0.861 

324.1 


TABLE II: The parameters and the saturation properties from MQMC-type models. The subscript 
M denotes a, u and p mesons for MQMC-S. For MQMC-MB, M refers to a; and p mesons. 
£>q ,/ 4 = 188.1 MeV and Zjv = 2.030 are used for all three cases. 


where 


ZT -M 


m 


M 


1 - 


a, 


EtLR. 


S(a) + 


E M 

^bag 


(25) 


The parameters that can produce the saturation at po are listed in Tab. |U| MQMC and 
MQMC-MB give us acceptable results for K and m* N , but the K value for MQMC-S is as 
large as twice of usually accepted ranges of values. To get a value of K within a reasonable 
range, we have tried many numerical searches over a wide range of parameter space, but no 
solution has been found so far. 


III. APPLICATION TO THE NEUTRON STAR 

We now apply the five nuclear models described in Sect. II to the neutron star mat¬ 
ter. Two basic assumptions of the neutron star matter are the charge neutrality and the 
(3— equilibrium. If we assume that only the nucleons and light leptons exist in the neutron 
star, charge neutrality is expressed as 

Pp = Y Ph (26) 

Z=e,/z 

where pi is the number density of particle i (= p, e, /i). Under (3— equilibrium, the processes 

n —> p + e~ + v e and p + e~ n + v e (27) 

occur at the same rate. The condition can be satisfied when the chemical potentials before 
and after the decay are the same. The chemical potential of each particle reads 

p n = \J + m% + guNUJ o - ^ g pN b 03 , ( 28 ) 
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(29) 

(30) 


lip = + mi) + S,i»i + -rUpNbm, 

li, = ^/kf + mf, 

where ki is the Fermi momentum of the lepton l (= e, /i). The chemical equilibrium condition 
is expressed as 

h-n = l~lp T ^ ) /P- (31) 

l 


A. EoS 


The EoS tells us the pressure (P) as a function of energy density (e), which are the 
diagonal elements of the energy-momentum tensor. The energy density and the pressure in 
the QHD model read 


£ = 


+ 



2 

03 


(32) 


P 


~\ m l ~ \ m N b (g aN v ) 3 - ( g aN a ) 4 + ^ml uj 2 + Kn 2 p 

1 ^ r kN k 4 1 f kl ki ,, 

—~ / / , =dk -\ -- } / , =dk. 

3?r N= n , p J v yjk 2 + m^ 2 3 tt 2 YP JW+rrf 




(33) 


The energy density and the pressure from the QHD-S model can be obtained from Eq. (TT?l) 
and Eq. El, respectively, by removing the cubic and quartic self-interaction terms of the cr- 
meson and replacing the free meson masses with the scaled ones given in Eq. ED). However 
the pressure from the QHD-S model has additional terms that stem from the density de¬ 
pendent meson masses and coupling constants. Thermodynamic consistency of these terms 
is discussed in Ref. [izj], and the explicit form is given as 
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FIG. 1: EoS from five different models. 


The energy density and the pressure of the MQMC-type models can be obtained by using 
QaN = 3 g q , g U N = 3 g% and g P N = g q p together with proper replacement of in-medium meson 
masses in Eqs. (13211 and (1341) . The resulting EoS curves for each model are shown in Fig.|l] If 
the pressure exerted by nuclear repulsion is strong, the matter becomes more incompressible, 
which corresponds to a large compression modulus and a stiff EoS. The EoS’s from the QHD, 
MQMC and MQMC-MB models whose IP s are close to each other (310.8, 285.5 and 324.1 
MeV, respectively) show similar behaviors at the energy densities considered here. 

The EoS from the scaling models (QHD-S and MQMC-S) turns out to be stiffer than the 
EoS from other models. This behavior can be understood by observing that the repulsion 
from the cn—meson is augmented at high densities (see Eq. © with replaced by m*) 
while the attraction caused by the a —meson is not strong enough to cancel the repulsion. 

To illustrate this argument clearly, we plot in Fig. [21 the magnitude of the contributions 
to the pressure P from N, a, uj and p. Each contribution P % (i = N, a, to, p) is defined as 


Pn 

Pa 


— y / , dk, 

N= n y* Jk 2 + m*N 

1 *2 — 2 T) I *2 — 2 t i 

2 m a <7 , Pu = -m u u> 0 , P p 



b 


2 

30 - 


The total pressure is approximately 


P — —P a + Poj + P p + Pn- 

The remaining terms such as the cubic and quartic self-interaction terms in QHD, extra 
terms for thermodynamic consistency and lepton contributions can be neglected since they 
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FIG. 2: Comparison of the magnitude of a , uj, p and the nucleon contributions to the pressure for 
the QHD (left) and the QHD-S (right). Definitions of Pi s are given in the text. 

do not determine the overall behavior of EoS. In the case of the QHD (see the left panel of 
Fig. I2J), Pa is non-negligible compared to P u in the density region considered here. Since 
cr-meson contributes to the pressure negatively, P a reduces the total pressure substantially, 
which in turn leads to a sizable softening of the EoS. On the other hand, P u from QHD-S 
is about 3 times larger than that of QHD, but P„ from QHD-S is more or less similar to 
that from the QHD. Thus the softening of the EoS due to P a in QHD-S is relatively weak as 
the density becomes high, and consequently the EoS from QHD-S is stiffer than that from 
QHD. 

Here, we need to remark that in Fig. [T] the EoS curves for MQMC-S and MQMC-MB 
terminate at a certain e. It is due to a breakdown of the SCO at a certain density. In 
Appendix, we show the details of why the SCC does not have a solution at some density. 

B. Composition 

The composition of the neutron star matter is represented by the number of particles 
divided by the total baryon number. The number of each species of particles is determined 
by the charge neutrality and the chemical equilibrium in the /3-decay. Muons can be created 
when the chemical equilibration condition between the electron and the muon, p e = p^ 
can be fulhlled. Given the values of a, u ) and &03 for each p, charge neutrality of Eq. (l26|) . 
chemical equilibrium of Eq. (EB, p e = p^ and p = p n + p p determine the number of particles 
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FIG. 3: The particle (to, p, e, p) fractions in neutron stars calculated by different nuclear models 
are plotted against neutron star matter densities p/po- Since the figure for MQMC-MB is similar 
to that for MQMC, we do not show it here. 


unambiguously. The results are shown in Fig. 01 

The compositions from the QHD and MQMC models are very similar to each other, but 
the number of the protons for the QHD-S or MQMC-S model is larger than that for the 
QHD or the MQMC, respectively. The ratio of the proton number to the neutron number, 
p p /p n is especially important in determining the cooling mechanism of the neutron star. 
Two major mechanisms of the cooling are emission of the neutrinos from the interior and 
radiation of the photons near the surface of the neutron star. The simplest neutrino emission 
mechanism is the so-called direct URCA (DU) process, which is nothing but the /3-decay as 
in Eq. (EH), and it is known as the most efficient neutrino-emitting process in the interior 
of the neutron star. In order for this process to occur, however, the ratio p p /p n should 
be larger than a critical value. The proton numbers from QHD-S and MQMC-S models 
increase faster than those from the QHD and MQMC models as the density increases. This 
implies that p p /p n from QHD-S or MQMC-S can reach the critical value at lower densities 
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and the /3-decay can take place over a wider region in the interior of the neutron star than 
the QHD or MQMC models predict, causing a rapid cooling of the neutron star. However, a 
recent work on the cooling of the neutron star |18] indicates that too low a threshold density 
for the DU process can result in an unrealistic cooling scenario. DU process can happen 
when both energy and momentum conservation is satisfied in Eq. m- Since neutrinos have 
energy of the thermal fluctuation around the Fermi momentum of the nucleon or electron, 
its momentum can be neglected. Then the momentum conservation condition for DU can 
be written as k n = k p + k e , where ki being the Fermi momentum of particle “i” . The 
threshold density pn at which DU starts to emit neutrinos is about po — 2.1p 0 f° r QHD and 
MQMC, and is about pr> — 1.6 p 0 for QHD-S and MQMC-S. These threshold densities are 


too low to explain the observation of the neutron star temperature 


Q. 


In a recent work 


the density dependent meson masses and coupling constants are considered within the 
framework of relativistic mean held theory. It is shown that they obtained lower value of pr, 
than that of our QHD model. Therefore the nuclear models we have considered in this work 
may need further refinement and adjustment of parameters to produce threshold densities 
high enough for the cooling process to be reasonable. 


C. Mass and radius of the neutron star 

The mass and the radius of the neutron star can be obtained by integrating Oppenheimer- 
Volkov equation with a given EoS. In Fig. 0 the neutron star mass with our EoS is plotted 
as a function of the central density e c (left panel) and of the radius (right panel) in units of 
solar mass M 0 . The solid line is the result from the QHD-S, the dashed line from the QHD, 
and the dotted line from the MQMC. QHD and MQMC give us very similar maximum mass 
M max ~ 2.0M 0 with R ~ 10.8 km and e c ~ 2.5 x 10 15 g/cm 3 . On the contrary, due to the 
stiff EoS the QHD-S gives us the maximum mass of about 2.9M 0 with the radius R ~ 13.7 
km and the central density £ c ~ 1.5 x 10 15 g/cm 3 . We could not extract the maximum mass 
of the neutron star for MQMC-S and MQMC-MB models since, as noted in the previous 
section, these models encounter a problem in solving the SCC at a density before reaching 
the maximum mass. 

Recent observations of the neutron star mass from the radio pulsars give us the mass 
range of the neutron stars as (0.8 2-2 )Mq R . The maximum mass of the QHD-S 
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FIG. 4: Neutron star mass as a function of the central density (left) and the radius (right). The 
solid line is the result from the QHD-S, the dashed line from the QHD, and the dotted line from 
the MQMC. 

is much larger than the upper limit of the observed values, so one may want to exclude 
the QHD-S among the models for the neutron stars. We may note, however, that more 
considerations are needed as follows. First, the QHD-S produces the saturation properties 
successfully. Secondly, we have considered only the nucleon degrees in this work, but other 
degrees of freedom can come in such as hyperons or the phases like meson condensation and 
quark matter. In general, when these exotic degrees of freedom appear, they reduce the 
Fermi momentum of nucleons through the /3-decay of the nucleons to exotic states. This 
results in a smaller contribution of Pn to the pressure, and thus the EoS becomes softer 
with exotic degrees of freedom than otherwise. Then the maximum mass of the neutron 
star becomes smaller. We have confirmed this behaviour by doing some more calculations 
including hyperons in the QHD and MQMC. The maximum mass of 1.7M 0 is obtained 
for QHD, and 1.5 M & for MQMC. It will be interesting to investigate how these maximum 
masses will change when medium effects on meson masses and coupling constants are taken 
into consideration. It is noted, however, that the inclusion of hyperons can give us too 
low a mass limit of the neutron star j^j], which is incompatible with the observations. In 
Ref. 0], a phase transition to quark matter is considered, and the maximum mass of the 
neutron star is obtained to lie in the range (1.4 ~ 1.7)M 0 . A transition to quark matter is 
one possibility, and at the same time transitions to other exotica are also possible. These 
possibilities should be considered for a better understanding. 
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IV. SUMMARY 


In this work, we have investigated the properties of the neutron star matter by including 
the effects of meson mass changes in medium. We have explored various models at the hadron 
(QHD and QHD-S) and at the quark (MQMC, MQMC-S and MQMC-MB) levels. We have 
incorporated in-medium meson masses in two ways: one by assuming a simple function 
of scaling in both QHD and MQMC models, and the other by treating heavy mesons as 
MIT bags within the framework of MQMC model. The EoS and the particle fractions for 
the neutron star matter are calculated. Scaling models show appreciable difference from 
the models with constant meson mass or meson bag models. The EoS’s from the scaling 
models are stiffer than those from other models. The scaling models show higher fractions of 
the protons than other models. Stiff EoS’s from the scaling models lead to heavy and large 
neutron stars whose maximum mass is about 2.9 M & which is larger than the observed values 
and the values predicted by other models by one solar mass. We may, however, note that 
the present work is just a first step toward the understanding of the relationship between 
the in-medium properties of physical quantities (mass, coupling constant and etc.) and the 
neutron star, and further studies need to be done. 
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Appendix 

In this section, we show a breakdown mechanism that occurs in solving the SCC of the 
QMC model. Since the SCC for MQMC-S and MQMC-MB are involved, we employ the 
QMC model for the purpose of a simple illustration. The SCC of the QMC model can be 
obtained by putting — 0 in Eq. ED, which leads to 

a = 3^- p s C N (a). (35) 

m l a 

The breakdown takes place when C'at(cx) changes the sign from positive to negative at a 
certain density. The numerator in the square bracket of Eq. ED causes the sign change. 
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A 

FIG. 5: The solid curve represents Eq. (El and the dots represent Q q = fi g (A) with Rm* obtained 
from the QMC calculation. In the left section divided by the solid curve, NU is greater than zero, 
which is allowed physically. QMC calculation has been done up to p = 5.3/>o- 


Rearranging the numerator of Eq. m, we obtain 

NU = n 2 q + 3Q 9 (Q 9 - 1 )Rm* - Z N Qfi, + (ft, - 1 )Rm*j (36) 

+ \( Rm q ) 2 + \' KRiB (^<? + “ l ) Rm i) > ( 37 ) 

where NU means the numerator. To show roughly how the sign of C^(a) changes, we 
make the following approximations. First, all the terms in Eq. El can be neglected since 
(Rm*) 2 and R A B are very small compared to the terms in Eq. (13611 in the density region 
we consider. Secondly, we use Zn = 2 instead of the actual value Zn = 2.030. With these 
approximations, we have 


NU ~ Qg(3 Rm* + 1) — Q q (5Rm* + 1) + 2 Rm* 


(3 Rm* + l)Q g 


2 Rm* 


(«, 


li¬ 


ps) 


Since (Q q — 1) is always positive, the sign change of Eq. (EHI) occurs at 

2 Rm* 

= - — 

3 Rm* + 1 

2A 

~ 3A - 1’ 


(39) 


where A = —Rm*. Since Q q is always greater than 1, the solution to Eq. EHD can exist for 
A > 1/3. By taking A as an independent variable, we plot in Fig. El Q q of Eq. (13H1) with 
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the solid line. Eq. © represents possible points at which the sign change occurs. The 
values of Q q obtained from the numerical solution of Eq. (1T5I) and the SCC of the QMC 
is shown in Fig. Q by the dots. In the actual numerical calculation with the QMC model 
(not MQMC model), of which we have not shown anything in this paper, solutions can be 
obtained up to p = 5.3p 0 at which a = 34.52 MeV, R = 0.574 and A = 0.533. In Fig. 0 
the two curves meet at around A = 0.527 (Q q = 1.812), which corresponds to the density 
p ~ 5.0po- If the small terms in Eq. © are included, the two curves will meet at around 
p = 5.3po, which is the breakdown point of the QMC in the actual numerical calculation. 
Similar breakdown behaviors are observed numerically for the MQMC-S and MQMC-MB 
models. For this reason, we could not calculate the EoS at higher densities, and consequently 
could not extract the mass and the radius of the neutron stars for these models. 
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